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Abstract — This paper studies the design and optimization of 
a limited feedback single-user system with multiple-antenna 
transmitter and single-antenna receiver. The design problem is 
cast in form of the minimizing the average transmission power 
at the base station subject to the user's outage probability con- 
straint. The optimization is over the user's channel quantization 
codebook and the transmission power control function at the base 
station. Our approach is based on fixing the outage scenarios 
in advance and transforming the design problem into a robust 
system design problem. We start by showing that uniformly 
quantizing the channel magnitude in dB scale is asymptotically 
optimal, regardless of the magnitude distribution function. We 
derive the optimal uniform (in dB) channel magnitude codebook 
and combine it with a spatially uniform channel direction 
codebook to arrive at a product channel quantization codebook. 
We then optimize such a product structure in the asymptotic 
regime of B — » oo, where B is the total number of quantization 
feedback bits. The paper shows that for channels in the real space, 
the asymptotically optimal number of direction quantization bits 
should be (M— 1)/2 times the number of magnitude quantization 
bits, where M is the number of base station antennas. We also 
show that the performance of the designed system approaches 
the performance of the perfect channel state information system 
as 2 _M + 1 . For complex channels, the number of magnitude and 
direction quantization bits are related by a factor of (M— 1) and 
the system performance scales as 2~ m as B — » oo. 

Index Terms — Beamforming, bit allocation, channel quanti- 
zation, limited feedback, multiple antennas, outage probability, 
power control. 

I. Introduction 

It is well established that the use of multiple antennas at base 
station can considerably improve the performance of wireless 
communication links. The realization of these improvements 
often requires channel state information (CSI) at the base 
station. In time division duplex (TDD) systems with uplink and 
downlink channel reciprocity, the base station can acquire this 
information via training on the uplink channels. In frequency 
division duplex (FDD) systems, however, the remote terminal 
needs to quantize and feedback the channel information to the 
base station. Such systems are generally referred to as limited 
feedback systems in the literature. 

The availability of CSI at the base station can significantly 
improve the performance of the limited feedback single- 
user and multiuser systems. In a single-user system, limited 
feedback of CSI provides the base station with power gain. For 
the multiuser multiplexing systems, the availability of CSI is 
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even more important. In such systems, the base station needs 
the users' channel state information to distinguish the users 
spatially and to perform rate and/or power control accordingly. 

A. Related Work 

1) Single-User Systems: Assuming perfect CSI, the ca- 
pacity of single-user channels with single antenna can be 
achieved by variable-rate coding at the base station, where 
the transmission rate is continuously adapted to the channel 
magnitude |TJ, Q- The practical implementation of variable- 
rate coding is achieved by adaptive modulation and coding 
(AMC), where the received signal-to-noise ratio (SNR) is 
compared with a certain set of thresholds and, based on this 
comparison, an appropriate modulation and coding index is 
chosen for the current channel realization |(3)-|(5). From the 
perspective of a limited feedback system design, these SNR 
thresholds are in essence the SNR quantization thresholds and 
effectively form a quantization codebook for the received SNR 
or equivalently for the channel magnitude. For an optimal 
system design, one therefore needs to optimize this codebook 
so that the average transmission rate is maximized for a given 
target decoding error probability at the receiver side JTJ, (2). 

Optimal design of limited feedback systems becomes more 
complicated when the stations are equipped with multiple an- 
tennas. The capacity of a limited feedback single-user channel 
with multiple antennas and temporally i.i.d. fading is shown to 
be achieved by sharing a codebook of transmission covariance 
matrices between the transmitter and receiver, where the 
receiver chooses the best transmit covariance matrix for the 
current channel realization and sends the corresponding index 
back to the transmitter ||6)-|[8|. Due to the complexities of the 
covariance codebook design, most of the work in this field 
focuses on either channel magnitude quantization j9)-fTT) or 
channel direction quantization fl2")-fT6), but not both. One 
can however easily show that for a limited-feedback system 
to approach the performance of the perfect CSI system, both 
channel direction and magnitude information are required fT7) , 
(TSJ. It is therefore essential to study the design of magnitude 
and direction codebooks jointly and derive the optimal split 
of feedback bits between the two codebooks. The problem 
of joint codebook design and optimization has not been well 
addressed and formulated in the literature and it is the main 
subject of this paper. 

2) Multiuser Systems: Although this paper is mainly con- 
cerned with single-user systems, it is informative to provide 
a short review of the literature on multiuser limited feedback 
systems. 

Single-antenna multiuser systems are widely addressed in 
the literature, e.g. (19), |20|. The problem in general is to 
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find a channel magnitude feedback mechanism that realizes 
the multiuser diversity gain, i.e. the scaling of the sum rate 
with the number of users in a large network of users. The 
design processes in these papers are however mainly numerical 
and they do not provide insight into the magnitude quan- 
tization codebook structure. The multiple-antenna multiuser 
systems are also widely studied in the literature [21 1-[26|. The 
availability of multiple antennas at the base station permits 
for spatial multiplexing of the users' data and for serving 
multiple users simultaneously. The work in (23) specifically 
shows that, for realizing the multiplexing gain, the number of 
feedback bits per user should scale logarithmically with SNR. 
In order to address the scheduling problem in large networks 
of users, a well-justified approach is to choose users with 
high channel magnitudes, low quantization errors, and almost 
orthogonal channel directions (22) , |26|. In order to realize 
the multiuser diversity gain in such networks, ]26| shows 
that one needs channel gain information (CGI) in addition 
to channel direction information (CDI); the CGI however is 
assumed to be unquantized in this work. The split of feedback 
bits between CGI and CDI quantization codebooks is studied 
in (27) . Although optimal bit allocation laws and codebook 
structures are not provided by [27 1 in a closed form, the 
authors are able to numerically show that as the number of 
users increase, more bits should be used for CGI quantization 
in order to benefit from multiuser diversity gain. 

B. Problem Formulation 

In this paper, we address the problem of channel magnitude 
and direction quantization codebook design and optimization 
for a multiple-input single-output (MISO) single-user channel. 

We are mainly interested in the optimal split of feedback 
bits between the two codebooks and also deriving the scaling 
of the system performance with the number of feedback bits. 
The system design problem is formulated as minimizing the 
average transmission power at the base station subject to the 
user's target outage probability. The optimization is over the 
user's channel quantization codebooks and the power control 
function at the base station. The problem formulation used 
here is typically more appropriate for fixed-rate delay-sensitive 
applications, e.g. voice over internet protocol (VoIP) and fixed- 
rate video streaming applications. The power control strategy 
that results from this formulation tries to fix the transmission 
rate by compensating for the channel fading and in this sense 
resembles the power control mechanisms used in Wideband 
Code Division Multiple Access (WCDMA) system standards 
[28 1 . Alternatively, it is also possible to fix downlink transmis- 
sion power and instead utilize adaptive modulation and coding 
to adapt the transmission rate to the user's channel quality as in 
Worldwide Interoperability for Microwave Access (WiMAX) 
and 3GPP Long Term Evolution (LTE) system standards |29|. 
Our problem formulation and results are not directly applicable 
to these latter models. The problem of optimal feedback bit 
split between channel magnitude and direction codebooks for 
such systems is not well studied in the literature and further 
analysis is required in this regard. 

The problem formulation used in our work is somewhat 



the power efficiency of limited feedback systems. The work 
in (32 1 specifically formulates the problem of minimizing the 
transmission power in a wireless sensor network subject to an 
average rate and average bit error rate constraint at the fusion 
center. The multiple single-antenna sensors in this work are 
assumed fully synchronized and therefore can be effectively 
modeled as a single multiple-antenna transmitter. Based on 
this model, the authors utilize adaptive modulation and coding 
and present a numerical approach to optimize the channel gain 
and direction quantization codebooks. The work however does 
not address the optimal split of the feedback bits between 
the codebooks and only numerically studies the scaling of the 
system performance with the number of feedback bits. 

Our work is also related to the literature on vector quan- 
tization codebook design, although the deign objective used 
here is quite different from the classical design objectives, 
e.g. minimizing the mean squared error (MSE) distortion 
measure. It should also be noted that the channel magnitude 
and direction codebooks used in our analysis effectively form 
a product quantization codebook for the user's channel vector. 
This product structure, also known as shape-gain quantization 



similar to the formulations used by |30|-[32| for improving 



1 33 1, provides several practical advantages including faster 
quantization and lower storage requirement for the quantiza- 
tion codebooks. In addition, most practical systems already 
have power control modules that essentially act as channel 
magnitude quantizers; therefore, a product codebook structure 
is more easily adopted in such systems. 

In order to study the problem of channel quantization 
codebook design, we first focus on the scalar channel mag- 
nitude quantization and characterize the optimal magnitude 
quantization codebook. It is shown that, as the codebook 
size N increases, the optimal quantization levels form an ap- 
proximately geometric sequence and hence become uniformly 
spaced in dB scale. Interestingly enough, the asymptotic 
optimality of uniform magnitude quantization (in dB scale) 
is not affected by the magnitude distribution function as long 
as certain regularity conditions are satisfied. This uniform in 
dB channel magnitude quantization codebook structure in fact 
provides a theoretical justification for using uniform (in dB) 
SNR quantization thresholds for power control in practical 
systems such as WCDMA and IEEE 802.1 In system standards 
©, (28). 

For the optimal uniform channel magnitude quantization 
codebooks, we are able to show that the gradient of the 
average transmission power (the objective function) with re- 
spect to uniform (in dB) quantization levels should diminish 
as 8(7V~ 3 / 2 ) as N — >• oo. We also derive a codebook that 
achieves such optimal scaling. 

We next form a product channel vector quantization code- 
book comprising the optimal uniform (in dB) magnitude 
quantization codebook and a spatially uniform direction quan- 
tization codebook and study its performance in the asymptotic 
regime of B — > oo, where B is total number of feedback 
quantization bits. After deriving the optimal power control 
and beamforming function at the base station, we optimize the 
product quantization codebook structure such that the average 
transmission power is minimized. We show that, for channels 
in real space, 
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1) The optimal number of direction quantization bits is 
M ~ 1 times the number of magnitude quantization bits, 
where M is the number of base station antennas. 

2) As B increases, the system performance approaches the 
perfect CSI system as 2~ W + I . 

The reason that we focus on real space channels in this 
paper, is to compare the real-space magnitude-direction bit 
allocation laws with similar laws in the quantization theory 
literature, which are based on conventional mean squared 
error (MSE) distortion measure. For channels in the complex 
space, one can apply the exact same approach to show that 
the number of magnitude and direction quantization bits are 
related by a factor of (M— 1) and that the system performance 
scales as 2~ m in the asymptotic regime of J3— »oo (see Section 
|IV-C] ). 

The remainder of this paper is organized as follows. Section 
[II] shows the asymptotic optimality of uniform (in dB) channel 
magnitude quantization. Section III presents the single-user 
system design problem in its general form and describes 
our approach in transforming it to a robust optimization 



problem. Section IV describes the product codebook structure 



and derives the asymptotically optimal magnitude-direction bit 
allocation law. Section [V] presents the numerical results and 
Section VI concludes the paper. 

Notations: All the computations in this paper are for the the 
real space. All logarithm functions throughput the paper are 
base 2. Also, the angle between any two unit-norm vectors 
u and v is defined as Z(u, v) = arccos|u T v| so that < 
Z(u,v) < vr/2. 

II. Channel Magnitude Quantization 

This section studies the structure of the channel magnitude 
quantization codebook. The resulting codebook structure is 
used later in the paper to form a product vector channel 
quantization codebook. As mentioned earlier, we assume real 
space channels throughout the paper. 

A. SISO Limited Feedback System Optimization 

Consider a single-input single-output (SISO) channel h £ 
R. The input signal s and the output signal r are related as 

r = ^/pjkjhs + z, (1) 

where E[|s| 2 ] = 1, P(h) is the transmission power, and z ~ 
Af(0, 1) is the Gaussian receiver noise. 

Assume a perfect CSI system, where the base station is 
required to guarantee a fixed target SNR j a at the receiver. The 
optimal power control strategy is therefore channel inversion, 
where the base station sets the transmission power P(h) 
according to 



P(h) 
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Now consider a limited feedback system with perfect CSI 
at the receiver and a feedback link from the receiver back to 
the base station with a capacity of B bits per fading block. We 
are interested in a system design that minimizes the average 
transmission power subject to a target outage probability q 
at the user side. The outage probability is defined as the 
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Fig. 1. The optimal transmission power P(h) is the step-like function 
shown in bold line. The channel magnitude quantization codebook Y = 

{y (1) ,j/ (2) ,--- ,y (Xr) } is Axed. 



probability that the received SNR drops below a target SNR 

To perform power control in a limited feedback system, 
the receiver quantizes the channel magnitude \h\ and sends 
the _B-bit index of the quantized magnitude level back to the 
transmitter; the transmitter then sets the transmission power 
according to the quantized magnitude. The limited feedback 
system design is therefore is a two-fold problem: 1) optimizing 
the power control function given the magnitude quantization 
codebook; 2) optimizing the magnitude quantization codebook 
itself. 

For the SISO system described above, the solution to 
the first problem is straightforward and can be expressed 
in a closed form as follows. Consider a channel magnitude 
quantization codebook 



Y 



{y w J 2) , 
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where N is the codebook size 1 . Here ?/ n )'s are the magnitude 
quantization levels that are used for quantizing the magnitude 
variable Y = \h\ 2 . For a fixed codebook Y of size N, we are 
interested in the power control function P(h) that minimizes 
the average transmission power for a fixed target outage 
probability q. With this objective, the optimal strategy is to 
accumulate all the outage scenarios in the lowest magnitude 
interval, i.e. the leftmost interval in Fig. [T] and set P(h) in 
other intervals to the minimum possible value that prevents 
outage. This minimizes the average transmission power since 
preventing outage with a higher channel magnitude requires 
less transmission power. 

Having the optimal transmission power in terms of the 
quantization codebook Y, we can now address the optimization 
of the quantization codebook itself. The objective is to find 
the codebook Y of a given size N that minimizes the average 
transmission power subject to a target outage probability q. 

Define the quantized magnitude Y as the lower end of the 
quantization interval that includes Y: 

f = y (n) if y {n) <Y <y {n+1 \ (3) 

where l<n<N and y(^ +1 )= f oo. The normalized average 

'The dot notation ( ) , throughput the paper, is used for all the variables 
associated with channel magnitude quantization. For the variables associated 
with channel direction quantization the double dot notation ( ) is used. 
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transmission power, according to Fig. [TJ is given by 



V(Y) = — E[P(h)} 
% 



E 



E 



(n) 



(4) 



where Q( n ) = F(y( n+1 )) -F (yW) and F(-) is the cumulative 
distribution function (cdf) of Y. 

According to Fig. Ill the first quantization level is fixed 
by the outage probability q: 



y 



(i) 



(5) 



where is the inverse cdf of Y. In order to optimize the 

codebook Y, we take the derivative of V(Y) with respect to 
y( n ) f or 2<n<N and set it equal to zero. This results in the 
following equation: 



f(y (n) ) 



l 



y(rt-l) y\ 

where /(•) is the probability density function (pdf) of Y. By 
approximating 

Q (n) * f(y {n W n+1) -y in) ), (7) 

for 2<n<N— 1, we achieve the following relation between the 
magnitude quantization levels: 



1 

/(«) 



)(«) 



(6) 



(n+l) 



y {n) 

An-l)- 



(8) 



i.e. the optimal quantization levels y( n \ 1 < n < N, form an 
approximate geometric sequence and therefore are uniformly 
spaced in dB scale. This result, surprisingly enough, does not 
depend on the distribution of Y = \h\ 2 . The geometric ratio of 
such a sequence however should depend both on the codebook 
size and the distribution of Y as it is shown next. 

B. Asymptotic Optimality of Uniform (in dB) Magnitude 
Quantization 

In order to rigorously show the efficiency of the geometric 
sequences for magnitude quantization, we proceed to bound 
the gradient of the objective function V in Q with respect to 
such quantization levels. By showing that the norm of the 
gradient goes to zero as A increases, we can then claim 
that such quantization codebooks are at least asymptotically 
optimal as N — > oo. 

Starting with some definitions, we define a geometric se- 
quence f' 9 '(r) with the ratio parameter r > 1 as follows: 



W = or™- 1 



1< n < N\ 



(9) 



where a = f y^ for notation convenience. We refer to these 
sequences as the uniform (in dB) magnitude quantization 
codebooks. In proving the asymptotic optimality of uniform 
(in dB) magnitude quantization codebooks, we will need 
the following regularity conditions on the channel magnitude 
distribution: 

Definition 1: Magnitude distribution regularity conditions: 
1) f(y) i s a positive differentiable function with bounded 
derivative over y > 0. 



2) r\ = lim^oo -f(y)/f(y) exists and t) ^ 0. 

3) E[Y] is bounded and therefore lim^oo yf(y) = 0. 
The logic behind these assumptions becomes clear in the 
discussions to come. 

In order to bound the gradient of the objective function 
V in Q, one can use the Taylor expansion to bound the 
approximation error in |7]) and obtain the following bound 
on the gradient: 

Lemma 1: For the uniform magnitude quantization code- 
book (r), we have 



\^Y>.!i)(r)'P\ 



1 

< 2 fl 



D, 



'N{r-lf 

where /i is an upper bound of f'(y) for all y > 0, and 

f(t)dt) 



D = 



\ ( i )///(//)- 



(10) 
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Proof: See Appendix [I] ■ 
As discussed earlier we want the gradient bound ( fT0) i to 
diminish to zero as the codebook size N increases. To achieve 
this, the geometric ratio r should be finely tuned with N, The 
idea is to set r appropriately, so that both terms on the right- 
hand side of (jTOf go to zero as N —> oo. Since r > 1, we can 
express r as a function of N as 

r = 1+C(N), 

where C(N) is an arbitrary positive function. Define 

C(A0 = -log£(A0/logA, 



so that 



Also, let 



1 + 



C(oo)= lim COO- 

JV->-oo 



in (10), both V N(r - l) 2 and D 
for D 



0, we need y . 



To achieve 1 1 VT 3 1 
should go to zero. From (TTT 
ar N-i qq_ Substituting r = 1 + A - ^^' and applying 
L'Hopital's rule, we get a necessary condition C(°°) ^ 1- 

Now, consider the first term in ( fT0| ), since C(°°) ^ 1 
and r = 1 + N~^( N \ we see that the first term cannot 
go to zero faster than 0(A~ 3 / 2 ). Since ||VP|| is bounded 
by the magnitude of the first term in (jTDJ, we conclude 
that in the class of geometric sequences, the best possible 
scaling for HVPH is 0(A -3 / 2 ), which is achieved only with 
r>l + e(A-!). 

C. Proposed Magnitude Quantization Codebook 

In the following, we propose a specific solution for r 
that achieves the optimal scaling. The uniform magnitude 
codebook corresponding to this value of r is referred to as 
the optimal uniform (in dB) codebook in the remainder of 
this paper. 

Definition 2: For any constant 0<c<oo, define the positive 
function C c (n) for n>l as the solution to the following 
equation: 



£ c (n) (1 + C c (n)Y 



(12) 
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Note that C c (n) is a well defined function, since the left-hand 
side of (jT2j is monotonic in C c (n) and has the range of (0, oo). 
Also define the function Cc( n ) — ~ log£ c (n)/logn, i.e. 

C c {n) = n-^ n \ 

It can be shown that for any 0<c<oo, 

lim ( c (n)=l. 



10' 



(13) 



(14) 



Definition 3: Define the optimal uniform magnitude quan- 
tization codebook Y* as 



Y* = Y ( ^(r*), 



(15) 



with 



r* = l + £ n/a (N), 

where a = F^ 1 {q) and t) = ]im y -+ 00 -f(y)/J'(y). Also 
define the corresponding quantized magnitude Y* similarly 
as in 0. 

Now the norm of the gradient with respect to Y* can be 
bounded as follows. 

Theorem 1: Assuming the regularity conditions, we have 
the following for the optimal uniform codebook Y* and any 

e>0: 



||V Y ^|| < ijitf-K/.W-*) +o(iV- 3 + e ) 



(16) 



where /i is defined in Lemma [T] 

Proof: See Appendix [j] ■ 
Now, considering the bound in ( fTo*) and noting that 
lirrin^oo ( c (n)—l, we have the following for any e > 0: 



|V Y *-P|| =0 N-i 



(17) 

which proves that the uniform codebook Y* is asymptotically 
optimal as N^oo and furthermore it obtains the best scaling 
possible within the class of uniform (in dB) codebooks. 
Such optimality is not affected by the channel magnitude 
distribution as long as the regularity conditions in Definition 
Q]are valid. 

It should be noted that the codebook that achieves the 
optimal scaling is not unique 2 . However, the function C C (N) 
in ( |T2] > with c = r\ja is specifically defined so that we achieve 
the residual scaling of o(N~ 3 ) in ( [16) . Moreover, this defini- 
tion provides a close approximation of the optimal magnitude 
quantization levels as shown in Fig. [2] The figure compares 
the uniform (in dB) quantization levels in Y* with the optimal 
quantization levels for a chi-square random variable Y with 
the pdf 

f ^-2¥rW) yf ~ le '^ (18) 

where the integer M>2 is the distribution parameter and T(-) 
is the gamma function. This is in fact the distribution of ||h|| 2 
for a M-antenna channel vector h = [hi, hi, ■ • ■ , /»jf] T where 
hi's are independent and distributed according to Af(0, 1). 
The optimal quantization levels are computed by numerical 
minimization of |4} with multiple random start points. 



"For example, the function C(N) 
achieve the scaling of 0(A r_3/ ' 2 ). 
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can also be shown to 
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Fig. 2. Uniform (in dB) magnitude quantization codebook ¥* vs. the optimal 
quantization levels for the chi-square random variable Y; Case 1) Af=3, 
q=10~ 5 , and iV=100; Case 2) M=4, (j=10- 4 , and Af=50 



In the remainder of this paper, we use the uniform codebook 
Y* for channel magnitude quantization. In order to study the 
performance of this codebook in terms of the corresponding 
average transmission power, we present an upper bound on the 
average normalized transmission power V(Y*) = E[l/Y*], 
where V(-) is defined in (ffl, and Y* is the quantized magni- 
tude for the optimal codebook Y* as described in Definition 

13 

First, assume a perfect CSI system where the transmission 
power P(h) is given by Q. The normalized average trans- 
mission power is therefore given by 



E 



" 1 " 




Y 


-j 

Jo 



1 

o V 



f(y)dy, 



(19) 



where the subscript SU stands for single-user. 

The following theorem presents an upper bound on 
V(Y*) =E[1/F*]. 

Theorem 2: Assuming the regularity conditions on the 
channel distribution function, we have the following for N>1: 



E 



where u> 



' 1 " 


< E 


' 1 " 






Y*. 




Y 



1 + N 



_ E[Y] 



(20) 



T ? ^SU.CSI 

Proof: See Appendix [I] ■ 
This result is used in the subsequent sections for optimizing 
the magnitude-direction product quantization codebooks. In 



Section III we present the general formulation for the optimal 
design of the limited feedback system and describe its robust 
variant. By considering the robust formulation, Section [IV] 
addresses the product channel quantization codebook opti- 
mization and presents the asymptotic bit allocation laws for 
the channel magnitude and direction quantization. 
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III. Vector Channel Quantization: The General 
Form 

Consider a single-user MISO channel heR M , where M is 
the number of base station antennas. The input signal s and 
the output signal r are related as 

r = y/P(h) h T v(h)s + z, (21) 

where P(h) and v(h) are the transmission power and the 
unit-norm transmit beamforming vector and z ~ AT(0, 1) is 
the receiver noise. The received SNR for this setup is given 
by P(h) |h T v(h)| 2 . 

Assume a limited feedback system with perfect CSI at the 
receiver and a feedback link with capacity B bits per fading 
block. The receiver quantizes the channel h and sends back the 
corresponding P-bit quantization region index. The transmitter 
then determines the transmission power and the beamforming 
vector using the quantized information. 

We are interested in a system design that minimizes the av- 
erage transmission power subject to a target outage probability 
q. The outage probability is defined as the probability that the 
received SNR drops below a target received SNR j g . 

Similar to the scalar channel quantization, the system design 
is a two-fold problem: 1) finding the optimal transmission 
power and beamforming vector given the channel quantization 
codebook; 2) optimization of the channel quantization code- 
book itself. In the following, we present the general form for 
such system design and present our approach in transforming 
it into a robust design problem. 

We start by some definitions: 

Definition 4: By a vector channel quantization codebook 
C of size N, we mean a partition of R M into N disjoint 
quantization regions S^ n \ l<n<N: 

C={5 (1) ,S (2) ,'",S m }. 

For every quantization codebook C, we also define a quanti- 
zation function 

5(h) : E M -> C, 

which returns the quantization region that hel* belongs 
to. 

For a given total number of feedback quantization bits B, 
the target SNR 7 , and the target outage probability q, the 
system design problem is formulated as follows: 



min E h [P(5(h))] 

P(S('h)), 
v(S(h)) 



s.t. 



prob 

\C\ = 



P(5(h))|h T v(5(h))| 
N = 2 B . 



<7o 



(22) 



<q, (23) 



where N is the quantization codebook size and the opti- 
mization is over the codebook C, the power control func- 
tion P(5(h)) : C — > R+, and the beamforming function 
v(5(h)) : C — > Hm, where Hm is the unit hypersphere in 
R M . 

An exact solution to this problem is intractable and the 



we fix an outage region of volume q and design the system 
such that the target SNR is guaranteed for all other regions, 
which we refer to as no-outage regions. With this approach, 
the outage probability constraint in ( p3| is replaced with SINR 
constraints that should be satisfied for all no-outage regions 
and therefore the problem in ( |22] i transforms to a robust 
optimization problem. 

To make this rigorous, define the outage region O C C such 
that prob [5(h) € O] = q and let the random variable 1(h) 
be the no-outage flag, i.e. 1(h) = l(S(h) G O c ), where the 
logic true function !(■) is 1 if its logical argument is true and 
it is otherwise. 

For a robust system design, we need to design the codebook, 
the power control function, and beamforming function in 
a robust manner such that the target SNR is guaranteed 
whenever /(h) = 1: 

min E h [P(5(h))] (24) 

P(S('h)), 
v(5(h)) 

s.t. inf P(5(h)) |w T v(5(h))| 2 > 7o J(h), 

for all het". (25) 

Note that by including the no-outage flag in the constraint 
( p5j ), this formulation returns P(5(h))=0 if the receiver is in 
outage, i.e. /(h) = 0. Also note that the outage constraint, 
prob [/(h) = 0] = q, is implicit in this formulation. 

It should be noted that the constraint ( f25| will not be feasible 
if the quantization region 5(h) includes the zero vector w = 
(unless /(h) = 0). This means that the outage region O should 
at least include the quantization region in C that encompasses 
the origin. 

For the robust design problem in |24|, the optimal power 
control function can be directly computed using the constraint 
as follows: 



P su d ^ f P(5(h)) 



V(h) 



inf |w^v(5(h))r 

wGS(h) 



(26) 



where the subscript SU stands for single-user. The problem in 
d24l) therefore simplifies to the following problem: 



min Eh 

C,v(5(h)) 



7o^(h) 



(27) 



numerical optimization approach is rather involved [17]. Our 
approach for solving this problem is a suboptimal one where 



_inf wes(h) |w T v(5(h))|' 

where the optimization is over the codebook C and the 
beamforming function v(5(h)). 

The problem in |27| cannot be solved without a specific 
assumption on the geometry of the quantization regions. In 
the next section, we further simplify this problem by assum- 
ing a magnitude-direction product structure for the channel 
quantization codebook. 

IV. Product Quantization Codebook Design and 
Optimization 

In this section we present the structure of the product chan- 
nel quantization codebook comprising a channel magnitude 
and a channel direction quantization codebook and study the 
optimal bit allocation between these two codebooks. 
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A. Product Codebook Structure 

1 ) Channel Magnitude Quantization: The structure of the 
channel magnitude quantization is exactly the same as the 
structure described in Section [TTJ with the difference that the 
magnitude variable is Y = ||h|| 2 instead of Y = \h\ 2 . 

Let C denote the set of quantization intervals for ||h|| = 
VY: 

C = {jWj (2) ,-,> ) }, (28) 



where =[y/y&) , \J yi n+x )). Here y^'s are the quantiza- 
tion levels for Y = ||h|| 2 , TV is the magnitude codebook size, 
and = f oo. Similar to Section Nil the first quantization 

level is fixed as follows: 



V 



W=F-\q), 



(29) 



where is the inverse cdf of Y = ||h|| 2 . 

According to Theorem [2j if we use the optimal uniform 
codebook Y* for magnitude quantization, we have the follow- 
ing for N>1: 



' 1 " 






^ Psu.csi ^ 



(l + + wtf-^/.W) , (30) 



E 



where p mcsi is defined in ( fl9] > with Y = ||h|| 2 . 

2) Channel Direction Quantization: We use a M- 
dimensional Grassmannian codebook U of size N for direction 
quantization: 



U 



{u«,u( 2 ),... ,u(*)}, 



(31) 



where u'™' vectors are M-dimensional unit-norm Grassman- 
nian codewords. 

The direction quantization regions are formed by mapping 
each channel vector h to a vector u(h) G U that has the 
smallest angle with h: 



u(h) = areminZ(h, u). 

uGO 



(32) 



The vector u(h) is referred to as the quantized direction for the 
channel realization h. The corresponding quantization regions, 



according to the Gilbert- Varshamov bound argument [34], can 
be covered by the following spherical caps: 



(33) 



where 



B (n) = {w e Um 



-(w,u 



< 



Here Hm is the unit hypersphere in H M , flW is the spherical 
cap around w- n > and 

(f> = arcsin S 

is the angular opening of the caps, where S is the minimum 
chordal distance of U. 

Covering the direction quantization regions with these 
spherical cap regions allows for a closed form solution for 
( p7| i. Moreover such enlargement of the regions increases the 
required transmission power and will result in an upper bound 
for the average transmission power. 



In order to describe the dependence between the angular 
opening <j> and the direction codebook size N, we need the 
following lemma: 

Lemma 2: For a M -dimensional real Grassmannian code- 
book of size N and minimum chordal distance 5, we have the 
following for large enough TV: 



5 < A\ M N J"- 1 , 



(34) 



where A» 



_ / y^r((M + l)/2) \ M-l 

^ M — \ r(M/2) J 
Proof: The proof is based on the Hamming bound as 
in fTJ) , with the difference that, unlike the complex space, 
closed-form expressions do not exist for the surface area of 
the real spherical caps (see Appendix [DJ. ■ 
We therefore have the following for large enough N: 



sin 4> < 4A M N 



(35) 



This bound, along with the bound in pO] ), is used in Section 
|I V-B | for optimizing the product codebook structure. 

3) Outage Region and Product Codebook Structure: Ac- 
cording to the robust design argument in Section [Til] we have 
to specify and fix an outage region O such that prob [h € O] = 
q, where q is the target outage probability. For this purpose, 
we set the outage region to be ball centered at origin: 



h llhl 



(36) 



where t/ 1 ) is the first quantization level, which is fixed 



according to 

Now, the product channel quantization regions are described 
as follows: 

C = (C x C) U O, (37) 



where C includes the magnitude quantization regions in ( |28j ) 
and C includes the direction quantization regions in ( [33) . 
Each quantization region in $7} , except the outage ball, is 
schematically a sector-type region as shown in Fig. [3] The 
variable Y in the figure is the quantized magnitude variable 
as defined in ([3]). Finally, noting ( [37] l, the size of the product 
quantization codebook C is related to the magnitude and 
direction quantization codebook sizes as follows: 



\C\ =NN+l. 



(38) 



Given the structure of the channel quantization codebook, 
the next section addresses the product codebook optimization 
and derives the bit allocation law between the channel mag- 
nitude and direction quantization codebooks. 

B. Product Codebook Optimization and the Optimal Bit Allo- 
cation Law 

As described in Section [Hi] we are interested in a robust 
system design that minimizes the average transmission power 
subject to a target outage probability as formulated in ( |24"| >. The 
system design includes optimizing the power control function, 
the beamforming function, and the quantization codebook 
structure itself. The optimal power control function is given 
by 



26b in Section III 



This section addresses the optimal 
beamforming function and the optimal product codebook 
structure. 
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Fig. 3. Sector-type product quantization region <S(h) for a channel realization 
h. 



First, fix the codebook C. Consider the optimization of 
the objective function in ( p7| ) over the beamforming function 
v(«S(h)). For a channel realization h G 5(h), the optimal 
transmit beamforming vector, according to (|27]>, is given by 



v(<S(h)) = arg max inf |w 2 

vGit M weS(h) 1 



(39) 



It can be shown that the optimal beamforming vector is in fact 
the quantized direction itself: 



v(<S(h)) - u(h). 



(40) 



This is expected noting the symmetry of the quantization 
region around u(h) as shown in Fig. [3] 

For the beamforming function v(«S(h)) = u(h), it is easy 
to verify that 



inf |w T v(S(h))| 



= inf |w T u(h)| = 7(h) cos 2 <b. 
wes(h) 1 1 

(41) 

By substituting this in d26| we have 



P. 



7qJ(h) 
f(h) cos 2 , 



(42) 



Note that 7(h)=l for h^C By taking the expectation of the 
both sides of |42} we obtain the following expression for the 
average transmission power: 



E[P SI 



7o 



-E 



(43) 



where Y is the quantized magnitude as defined in ([3]). 

We are now ready to optimize the product codebook for 
robust system design. For this purpose, we present an upper 
bound for E[P SU ] in ( |43] l and use it to formulate the product 
codebook optimization. The computations are asymptotic in 
the codebook sizes, i.e. we assume N, N 3> 1 throughout the 
optimization 

The expression E[l/Y] in ( |4"3"j ) can be minimized by using 
the optimal uniform (in dB) magnitude codebook Y*. By 
setting Y = Y* in ( |43| ) and using the bound in ([30} and 
also the bound on <f> in (|3~5j), we obtain the following upper 



bound for the average transmission power: 



p 



< 



l + JV-Cn/.(*) +w JV-2C,/.(*) 



(a) 



1 

1 + N~ 



16\ 2 M N-^r 



1 - 16X 2 M N-v^ 



(b) 



where P„ 



def 



according to ( fl4| ) 



N^ 1 + 16\ 2 m N-t^, (44) 
SI . The approximation (a) holds since 



lim C v/a (N) 



1. 



Also the approximation (b) holds since we assume N,N 1. 

For a total number of feedback bits B, we have the 
following constraint on the magnitude and direction codebook 
sizes: 

NN + 1 = N = 2 B , 

where N is the size of the product quantization codebook 
C. Since the computations are asymptotic in N and N, we 
approximate this constraint as 

NN = 2 B . 

The codebook size optimization therefore simplifies to the 
following problem: 



min N^ 1 + 16\ 2 M N~ t^t 

N,N 

s.t. NN = 2 



(45) 



B 



The following theorem gives the optimal channel magnitude 
and direction codebook sizes: 

* def ' " def 

Theorem 3: Let B = \og{N) and B = log (TV) denote re- 
spectively the number of quantization bits assigned to channel 
magnitude and channel direction quantization. By solving ( |45| ) 
the optimal number of quantization bits are given by 



B= 



-B — k„ 



M + l 
•• M - 1 



(46) 



(47) 



where k - M=l ] n(r 32A m 
wnere k su — M+1 log M1 . 

Proof: The optimal number of bits are derived simply by 

applying the Lagrange multipliers method to the problem in 

(g3]>. ■ 

Corollary 1: For a single user system with M antennas 

at the base station, the number of magnitude and direction 

quantization bits are related as follows: 



•• M — 1 ■ M + l 

B = — - — B H -r— K c 



(48) 



2 ~ ' 2 su ' 
where k su is defined in Theorem [3] In the asymptotic regime 
of B — > oo, the optimal number of magnitude and direction 
quantization bits are therefore related as 



•■ M - 1 • 

B = B. 



(49) 



It should be noted that the codebook optimization problem, 
as described in this section, can be viewed as a classical 
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constrained quantization codebook design problem if one 
considers the average transmission power in (j24j» as the quan- 
tization distortion measure. The relative bit allocation law in 
( |49| i can specifically compared with the results in p5| , which 
studies the optimal shape (channel direction) and gain (channel 
magnitude) quantization with the mean squared error (MSE) 
distortion measure. The authors show that, with the MSE 
measure, the asymptotic number of shape-gain quantization 
bits are related as 

B = (M - 1)B. 

Comparing this with ( |49| l implies that the number of direction 
quantization bits for a product codebook optimized for the 
MSE measure differs by a factor of two. 

Theorem [3] is the final part of the product codebook 
optimization for the single-user system. By combining this 
structure with the beamforming and power control functions 
in ( |40] > and ( |42) , we achieve a robust system design for a 
limited-feedback single-user MISO system. 

Finally, in order to derive the performance scaling of the 
designed system with the feedback rate B, we substitute the 
optimal codebook sizes of Theorem [3] into ( |44| i and compute 
the average transmission power. 

Theorem 4: For a single-user system with M antennas at 
the base station and a total feedback rate of B bits per fading 
block, we have the following asymptotic upper bound as B — >• 
oo: 

nP sv ] < P wm (l + , (50) 

where P SUCSI = 7 P S u,cs, and 

_ 16(M + 1) / ygM - l)r((Af + l)/2) \ 2/(M+1) 
<7su ~ M - 1 V 32r(M/2) J 

Proof: The proof uses basic calculations, which are 
omitted due to the space limits. ■ 
Thus, as the feedback rate B increases, the performance of 
the designed system approaches the performance of the perfect 
CSI system as 2~*?+t. 

C. Complex Space Channels 

As mentioned earlier, the results in Section |IV-B| are based 
on the assumption of real space channels. The same approach 
however can be used to derive the bit allocation laws for 
complex space channels. The only main difference is in the use 
of Hamming bound in p4) i. For the complex space channels, 
the Hamming bound is as follows JT3): 

5<2N^^^. (51) 

By using this bound instead of |34]i, the remainder of the 
analysis can be applied in a similar fashion to derive the 
optimal magnitude-direction bit allocation law. By doing so, 
one can show that for complex channels and in the asymptotic 
regime of B — >• oo, the number of magnitude and direction 
quantization bits are related as B={M—1)B instead of the law 
in (|49j for real channels fTH) . Furthermore, for the complex 
space, the the system performance scales with the feedback 

_ B_ _ 2B 

rate as 2 a/ instead of 2 M + 1 for the real space. 




10° 10 1 10 2 10 3 

Direction quantization codebook size 



Fig. 4. Sine of angular opening <j> as a function of direction codebook size 
TV for Grassmannian codebooks in R 5 . 

V. Numerical Results 

In the process of deriving optimal bit allocation laws in 
Theorem [3] we used upper bounds in ( fl6| i and p5] l and 
the approximations in ( |44] > to find an analytically tractable 
approximation of the original objective function in ( |43"] >. In 
order to examine the accuracy of the magnitude and direction 
bit allocation laws in Theorem [3] this section compares these 
closed-form solutions with the corresponding numerically op- 
timized bit allocations. 

In order to find the optimal bit allocations, we need to 
numerically compute and minimize the average transmission 
power in ( |43]l in terms of the magnitude and direction code- 
book sizes N and N. The minimization is subject to the 
constraint NN=2 B or equivalently B+B=B, where B and 
B are the number of bits assigned to magnitude and direction 
quantization and B is the total number of feedback bits. 

The expression for the average transmission power in ( |43"j ) 
comprises two terms that are controlled independently by the 
magnitude and direction quantization codebooks. The term 
E[l/y] in ( |43| l is controlled by the channel magnitude quan- 
tizer and depends on the magnitude quantization codebook 
size N. The variable on the other hand is controlled by the 
direction quantizer and depends on the direction codebook size 
N. 

For the channel magnitude quantization, we use the uniform 
(in dB) codebook Y* defined in Definition [3] and numerically 
compute the term E[l/Y] using Q. In our numerical results, 
we use the chi-square distribution in ( fTS) ) with M = 5 as the 
channel magnitude distribution. For the direction quantization, 
on the other hand, we use real-space Grassmannian codebooks 
and rely on the available numerical tables that give the 
angular opening as a function of the codebook size N |36|. 
The values of sin</> for channel space dimension M=5 and 
codebook sizes up to iV=100 are shown in Fig. [4] 

Since the numerical tables for Grassmannian codebooks are 
only available for moderate values of N, we need to derive 
an extrapolation of the available data for larger values of N. 



to 



20 
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Total number of feedback bits 



Fig. 5. Magnitude and direction bit allocations for M = 5 antennas and 
target outage probability q = 10 — 4 . 
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Total number of feedback bits 

Fig. 6. Distortion measure T>(B) for M = 5 antennas and target outage 
probability q = 10 -4 . 



For this purpose we fit a line in the least-squares sense to the 
available data (N, sin tfi) in the logarithmic scale. Fig. [4] shows 
a line fitting that uses the data points with 80 < N < 100 as 
input. The slope of this line is forced to be — M 1 _ 1 in order to 
match the Hamming upper bound in ([35]> 3 . This leads to the 



following approximation of sm<fi for codebook sizes N > 100: 

sin^w 1.9833JV~^, (52) 

where M = 5 is the channel space dimension or equivalently 
the number of base station antennas. In numerical computation 
of the objective function in ( |43j ), we rely on the approximation 
in ((52]) whenever N > 100. 

To find the optimal bit allocation for a given number of 
feedback bits B, we numerically compute and compare the 
objective function |43) for different integer pairs (B,B) that 
satisfy B+B=B and choose the best pair. The resulting 
bit allocation is shown in Fig. [5] where the target outage 
probability is set to q = 10~ 4 . The figure also compares the 
results with the values of B and B in ( ftp*) and ( |47| ), which are 
rounded to the closest integer numbers. As the figure shows 
the analytically optimized bit allocations in Theorem [3] are 
at most two bits away from the numerically optimized ones. 
Furthermore, as the number of feedback bits B increases, the 



asymptotic result of B 



M-i 



B becomes more accurate. 



In order the evaluate the system performance under these 
bit allocations, we define a distortion measure as follows: 



V{B) 



p 



(53) 



where E[P SU ] is the average transmission power with limited 

= 7oP S u,cs, is the average 



feedback as defined in (R31 and P. 



3 According to the Gilbert- Varshamov bound |l3| , (34J, we have 

sin^aiV- 1 /^- 1 ), 

for some constant a > 0. By combining this bound with the Hamming bound 
in (35) , it is evident that s'm(f> should asymptotically scale as N~ l) 
with the codebook size N. The slope of the fitted line in logarithmic scale 
should therefore be equal to — m—1 ' 



transmission power with perfect CSI as defined in Theorem |4] 
and 7 is the target SNR. For the channel magnitude with chi- 
square distribution in (|18|l, one can show that p svc , 



M-2 

and therefore P SUCSI = j^. 

The distortion measure T>{B) clearly depends on the total 
number of feedback bits B and also on how these bits are 
split between the channel magnitude and direction quantization 
codebooks. It should be noted however that the distortion mea- 
sure does not depend on the target SNR, since 7 Q appears as a 
multiplicative factor both in the enumerator and denominator. 
As a result, the optimal magnitude-direction bit allocation, 
which minimizes the distortion or equivalently E[P SU ], does 
not depend on the target SNR as it is verified by Theorem [3] 

Fig. [6] compares the distortion measure when one uses the 
numerically optimized bit allocations versus the case where 
bit allocations in |46} and ( |47| ) are used. As it is shown in 
the figure, the two bit allocations show a close performance, 
which verifies the accuracy of the closed-form bit allocation 
laws in Theorem [3] Fig. [6] also shows the distortion measure 
applied to the upper bound in Theorem [4] which leads to the 
following asymptotic upper bound: 



V{B) < a su T 



(54) 



Although the quantization codebook design objective used 
in our work, i.e. minimization of the average transmission 
power subject to a target outage probability, is quite different 
from those used in standard quantization theory, the distortion 
measure upper bound in |54} still resembles the asymptotic 
distortion bounds derived in the literature of high resolution 
quantization theory |33) , (37j. The authors of |37[ specifically 
develop a rather general framework for analysis of high 
resolution vector quantization with locally quadratic distortion 
functions and show that the average distortion function scales 
asymptotically with the number of quantization bits B as 
2"t, where k is the degrees of freedom of the vector under 
quantization. One would therefore naturally expect a denom- 
inator of M instead of M+l in the exponent of the upper 
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bound in ( |54| ), since the channel vector is M -dimensional. 
This difference, we believe, is due to the difference between 
the distortion measure used in this paper with the standard 
quadratic (e.g. MSE) distortion measures used in the literature. 
Nonetheless, it is still possible to present an intuitive justifi- 
cation, although inevitably inaccurate one, for the difference 
in the asymptotic scaling. We describe this justification next. 

We first note that for a M-dimensional channel vector 
in real space, the channel magnitude is a one-dimensional 
variable. The channel direction, on the other hand, is M— 1 
dimensional due to its unit-norm constraint. Now, according 
to the asymptotic bit allocation law in ( |49| ), every bit assigned 
to channel magnitude quantization translates to two bits for 
channel direction quantization. It is therefore possible to say 
that in the process of dividing the available number of bits 
B among the available degrees of freedom, the degrees of 
freedom in the channel magnitude space counts twice the 
degrees of freedom in the channel direction space. Hence, we 
can assume an equivalent degrees of freedom and calculate it 
as (M - 1) + 2 = M + 1. 

A similar interpretation can also be applied to channel 
quantization in complex space. In this case, the channel is 
2M dimensional, the channel magnitude is one-dimensional, 
and the channel direction is 2M— 1 dimensional. Nevertheless, 
as far as the power control and beamforming problem in this 
paper is concerned, any channel vector h is equivalent to all 
the channels in the form of e j9 h, where 6 is an arbitrary 
real number. This means that the channel direction space has 
effectively 2M — 2 degrees of freedom. Now by following the 
intuition described for real space channels, we can assume an 
equivalent degrees of freedom for channel quantization and 
calculate it as (2M - 2) + 2 = 2M. This would lead to an 
asymptotic scaling of 2~™ — 2~m, which coincides with 



the scaling mentioned in Section IV-C for complex channels 



VI. Conclusions 

This paper studies the design of single-user MISO system 
with limited CSI at the base station. The problem is formulated 
as the minimization of the average base station transmission 
power subject to the outage probability constraint at the 
user side. We show that the asymptotically optimal channel 
magnitude quantization codebook is uniform in dB scale. This 
result does not depend on the channel magnitude distribution 
function as long as some regularity conditions are satisfied. 
Combining the uniform in dB magnitude codebook with a 
spatially uniform channel direction quantization codebook, we 
form a product channel quantization codebook and optimize 
the quantization bit allocation as B — » oo, where B is the 
total number of feedback quantization bits. It is shown that 
for channels in real space, as B — > oo, the optimal number 
of direction quantization bits is M ~ 1 times the number of 
magnitude quantization bits, where M is the number of BS 
antennas. We also show that, as B increases, the performance 
of the designed system approaches the performance of perfect 
CSI system as 2 M + 1 . For channels in complex space, on the 
other hand, the number of magnitude and direction quantiza- 
tion bits are related by a factor of (A/— 1) and the system 
performance scales as 2~*r as B — > oo. 
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Appendix I 

For notation convenience, the notations y^ n > and QW in 
Section [IT] are replaced with y n and Q n . 



A. Proof of Lemma [7] 

Proof: We first note that y 1 =a=F~ 1 (q) is fixed accord- 
ing to |5|. By taking the derivative of Q with respect to y n 
for 2<n<N, we have 



dV 

dy n 



f{Vn) 



1 



1 



Vn-l Vn 



Q n 



(55) 



For 2<n<N, we have the following by adding and sub- 
tracting a fixed term to (55) : 



dV 

dy n 



=f(y n ) 

+ f(Vn) 



1 



1 



Vn-l Vn 
Un+l ~ Vn 
Vn 



f(Vn) 
Qn 

vl ' 



y n +i - y n 



in, 



(56) 



The first two terms in ( |56) cancel out since y n = ar" _1 
form a geometric sequence. By using the definition of Q n , 
we therefore have the following for 2<n<N — 1: 



dV _ 1 

dyn vl 



f(Vn) {Vn+1 ~ Vn) 



f{y) dy 



(57) 



By applying Taylor's expansion to the function g(t) 
So f(y) dy for y n <t<y n+1 , we have 



f(y) dy - f(y„) (y„+-L - y n ) 



Vn 



(58) 

where /x is an upper bound on g"(t) = f'(t) for y ra <i<y n +i- 
By combining ( fST) and ( |58] l, we therefore have the follow- 
ing for 2<n<N - 1: 



dV 



< \»{y n+ i - Vnflvl = \rir - l) 2 - (59) 



Now for the last quantization level y n , n = N, we have 
dV 



dy n 



f{Vn) 

1 

y; 
i 

v 



i 



Vn-l Vn 



VnfiVn) 



1 " 




On 


Vn_ 




vl 


Vr 




- 1 


Vn- 


-1 





l - ({r-l)y n f{y n )- j™ f{t)&tj 



=D. (60) 



n=N 
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By combining d59} and (|60|) we therefore have: 



I 

Y 



' N-l 



_''i(r)^ , ii= ( 

n=2 
'jV-1 

E 



< 



n=2 



dV 



&P 



9y 



N 



D 























l 






A 



Y 



<-nVN {r-lf + D. 



(61) 



B. Proof of Theorem [7J 

Proof: According to (15) , the geometric sequence param 
eter for Y* is given by 

r* = l + £„ /a (#) = 



(62) 



where 



y^oo f'(y) 



By substituting the value of r* in the upper bound of Lemma 
[TJin ( |6T) , we have 

||Vy*7>|| < Li/H 2 ^*)-*) + D, (63) 



where 



D = 



(r-l)y„/(y„)- / /(t)di 



(64) 



n—N 



Noting the definition of r* in (|62j), the definition of the 
function C c (n) in ( fT2) , and also that = ar^ -1 , we have 



(r* - l) yxr - £, /a (ff) a (l + £,/#)) 
therefore 

^-^L. =„#&,/.(*). 
r* — 1 



0-77/0 = 77, 
(65) 

(66) 



Note that according to ( |T4] >, li m jv-j.oo C^/aC^) = 1» an d 
therefore y^ —> 00 as N —> 00. 

In order to bound D in d63l, we substitute d65|l in d64l): 



D-- 



I{Vn) 

' y 2 



, yjv/(^jv) 

4 



1 



/(yjv) 



/(*)dt 



1 



f(t)dt 



N) Jy$ 



(67) 



Now noting that y^ — > 00 as N — > 00, we have 



JV 



Inn [°° /(t)dt= lim -p^ = 7?. (68) 



/(y 



Moreover, according third regularity condition in Definition [JJ 
in Section ITT1 we have 



lim yfrfiyfr) = 0. 

JV->-oo 



Fig. 7. Proof of Theorem [5] 



Combining ( |68[ i and ( |69| l with ( |67| i, we have the following 

as N — > 00: 



D = (l/4)=o(iV- 3 ^)). 



(70) 



And since lim^^^ ( v / a {N) = 1, we have Z? = o(iV~ 3+e ), 
for any arbitrary e > 0. By substituting this in ( |63~| , the proof 
is complete. 



C. Proof of Theorem [2] 
fVoq/: 

Consider the shaded area A in Fig. [7] We clearly have 

^su.csi ■ 

(71) 

On the other hand, according to Fig. |7j 

A= y = ly = 1 ( v( y*) - 9i 

By combining ( f/T) and ( |72") i, we have 

7>(Y*) < r> su , CSI + 9± = (1 + /V-^/.W) ftuca + 
y at > ' 



r 1 1 

.A < / -/(y) dy = E[-] = ftl 
Jo J * 



(72) 



Or 
Vn 

where we have used the definition of r* in ((62 
Now from Markov's inequality. 



= prob[F > y^] < 



and therefore by using d66|, we have 



E[Y] 
Vn ' 



QN < m_W]<r-2C v/a (N) 
By combining f73] > and f75) >, we have 



(74) 



(75) 



(69) where u) 



Psu,csi 

_ E[y] 



VPs 
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Appendix II 
Proof of Lemma[2] 

Let iljvf denote the unit hypersphere in R M . Also let U 
denote a Grassmannian codebook of size N and define 5 as 
the minimum chordal distance of this codebook, as defined 
in fl3] . Consider an arbitrary unit vector u and define an 
spherical cap around it as 



where 



Bt/> = {w £ tt M |Z(w,u) <tp} 



ip — arcsin 5/2. 



By applying the Hamming -bound argument, as in fl3) , an 
upper bound on the Grassmannian codebook size can be 
achieved as follows: 



N < 



A(il 



M ) 



A(^) ' 



where A(-) denotes the surface area of its argument. 
The area of the unit hypersphere in R M is given by 



A(lt 



MS 



NIC 



M ■ 



(76) 



(77) 



where Cm = r(M/2+i) ' an( ^ ^(') * s tne g amma function. The 
area of the spherical cap on the other hand is given by 



A(B^,) = 2{M — 1)Cm-i [ sm 

Jo 



M-2 



<p dip. (78) 



Since the minimum chordal distance of the codebook 
tends to zero as codebook size tends to infinity, we have 
lim^^^ ip = 0. Therefore, for < ip < ip and any < e < 1, 
we have the following for large enough N: 



,M-2 



„M-2 



i.M-1 



: tp> (l-e)(/r 
By using this inequality in ( f78] l, we have 
A(B^) > 2(1 - e)C M -iV 
By combining f76] i, ( |77] i, and ( f80] >, we have 

^ 2(1 - e)C M -Vp M - 1 ' 
and therefore by using ^ = arcsin 5/2, we achieve 

((l-eW^-^a-eXV^-Ml 



(79) 



(80) 



(81) 



e) sin* 7 " 1 V 



<(l-ety 



M-l 



< 



MC 



M AT-l 



2Cm-i 



AT 



(82) 



We therefore have the following for large enough N and any 

< e < 1: 

2 •• i 
5< \ m N~t^, (83) 



1 - e 



where 



/ MC 



V2C. 



M-l 



/iT((M + l)/2) 

r(M/2) 



By setting e = | in ( |83"j ), we achieve the bound in Lemma [2] 
and the proof is complete. 
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